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Abstract 

A precise analysis of the pp elastic scattering data at 7 TeV in terms of its amplitudes 
is performed as an extension of previous studies for lower energies. Slopes Br and Bj of 
the real and imaginary amplitudes are independent quantities, and a proper expression 
for the Coulomb phase is used. The real and imaginary amplitudes are fully disentangled, 
consistently with forward dispersion relations for amplitudes and for slopes. We present 
analytic expressions for the amplitudes that cover all t range completely, while values of 
total cross section a, ratio p, Bj and Br enter consistently to describe forward scattering. 
It is stressed that the identification of the amplitudes is an essential step for the description 
of elastic scattering, and pointed out the importance of the experimental investigation of 
the transition range from non-perturbative to perturbative dynamics, that may confirm 
the three gluon exchange mechanism observed at lower energies. 
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I. INTRODUCTION 



The determination of the imaginary and real parts of the amplitudes of elastic pp and 
pp scattering is an essential step for the understanding of the underlying dynamics 
governing hadronic interactions and for the construction and critical evaluation of 
theoretical models. This question has now increased renewed interest in view of 
the experimental activity at new energy frontiers of the LHC, particularly with the 
forward physics measurements of the TOTEM Collaboration [U |2], and also with 
the possibility of measurements in the large \t\ range. However, the identification of 
amplitudes is not at all a trivial task. There is no methodology completely free from 
subjective judgment and of interpolations and extrapolations based on investigation 
of data at all energies. 

Efforts to identify amplitudes were made in a treatment originated in the Stochastic 
Vacuum Model [3], resulting in good description of the data on da/dt at all energies 
from 20 to 1800 GeV [I]. In particular, it was made clear to establish the regularity 
in the general features of the real and imaginary amplitudes, with identification of 
magnitudes, signs and zeros of the individual parts. 

The constraints imposed by the well-known dispersion relations for the amplitudes 
and also by the new dispersion relations for the slopes [5] help to control the forward 
scattering parameters, and are an essential tool for the construction of valid analytic 
representations of the data. In the present work, we extend the previous analysis to 
the recent pp elastic scattering measurements of the LHC TOTEM Collaboration 
[U |2] at y/s = 7 TeV, obtaining a very precise description of the amplitudes and 
of the differential cross section in the whole \t\ range, including location, depth and 
shape of dip and bump. The quantities a , p and slopes Bj and Br that are the 
common parameters of the basic features of the forward scattering data, are also 
determined precisely from this all-t analysis. Furthermore, we show that the 7 TeV 
data favor the assumption of universality of the contribution from the perturbative 
three gluon exchange process at large \t\ that has been observed at 27.4 GeV [6] and 
confirmed at higher ISR energies. 

Our analysis shows why the dip/bump or inflection structure, that appears in pp 
and pp elastic differential cross section at ISR and Tevatron energies, has also been 
observed at 7 TeV [U |2] . This is due to the monotonic displacement of the zero of the 
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imaginary part as the incident energy increases, approaching the first zero of the real 
amplitude, with formation of a marked dip at about 0.5 GeV 2 . We remark that at the 
energies 540 and 1800 GeV no dip-bump appears because the imaginary zero is in the 
middle between the two zeros of the real part, with compensating contributions. This 
analysis also suggests that, mainly as consequence of the continued displacement of 
the zero of the imaginary part, an even more pronounced dip/bump structure should 
appear for 14 TeV pp scattering at a slightly smaller value \t\ ~ 0.4 GeV 2 . 

We organize the present paper as follows. In Sec. [XXJ we introduce the representation 
of our amplitudes. This representation contains the exponential nature of the very 
forward domain, with associated Coulomb amplitude and phase, and the larger \t\ 
behavior with dip and bump, together with the contribution from the perturbative 
three-gluon exchange process, all in unique analytic forms. We also show the form 
of the nuclear amplitude in impact parameter space, originated in studies with 
the Stochastic Vacuum Model. In Sec. |III| we show our results for the scattering 
amplitudes obtained from a precise fit to the experimental data at 7 TeV. It is shown 
that the forward scattering parameters determined from our analytic expression for 
all t values are in agreement with those determined by the TOTEM collaboration [2] 



. In Sec. [TV] we discuss the physical nature of the amplitudes, the position of zeros 
and its relation to the observed dip/bump structure. These properties are compared 
to the 52.8 GeV case, where a wider t domain has been measured. To make clear the 
physical meaning of the properties of the amplitudes, a comparison with another 
model is made. Finally in Sec. [V] we give a summary of the present work, and 
discuss the consequences of our analysis, together with the future perspectives. In 
Appendix, we derive the Coulomb phase function extended to deal with the case of 
different Bj and B R slopes . 



II. ANALYTIC REPRESENTATION OF AMPLITUDES OF pp ELASTIC 
SCATTERING 

We write the differential cross section as 

^ t =(hcf \T R (s,t)+ l T I (s,t)\ 2 . (2.1) 
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In the following discussion, we use the unit system where a is in mb (milibarns) and 
the energy in GeV, so that {he) 2 = 0.3894 mb GeV 2 . 

Tr and Tj, with dimensions GeV -1 , contain the nuclear and the Coulomb parts in 
the forms 

T R (s, t) = (s, t) + V^F c {t) cos(«$) , (2.2) 

and 

T/(s, t) = Tj N (s, t) + VnF c \t) sin(a$) . (2.3) 

The Coulomb scattering amplitude F c '(s,t) enters in the pp/pp amplitudes with 
the form 

F c (t) <**M = (_/+) F 2 roton (t) , (2.4) 

where a is the fine-structure constant, $(s,t) is the Coulomb phase and the proton 
form factor is written 

FprotonW = [0.71/(0.71 + \t\)} 2 . (2.5) 
The total cross section is given by 

a = 4^ (he) 2 Tj N (s } t = 0) . (2.6) 

Tft(s,t) and T/ V (s,t), respectively the real and imaginary parts of the properly 
normalized scattering amplitude of the strong interaction, are smooth and regular 
functions of s and t, and the Coulomb amplitude is relevant in the very forward 
range \t\ < 10~ 2 GeV 2 . We neglect spin effects. 

A. Nuclear Amplitudes for all-t Values 

To obtain precise description of the elastic da/dt data for all \t\, we use the forms 
that have been introduced before and shown to be successful with ISR and Fermilab 
data |4J, together with the assumption of the perturbative three-gluon exchange 
amplitude |7j, writing 

Tg(s,t) = a K (s)e- Ms)W + X K (s)^ K ( lK (s),t) + 5 K , R R ggg (t) , (2.7) 

where ^k(jk(s), t) are shape functions, described below, and R ggg (t) represents the 
contribution from the perturbative three-gluon exchange amplitude. The label K 
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means either K = R for the real amplitude or K = I for the imaginary amplitude, 
and 8k,r is the Kronecker's delta symbol, that is, the last term only contributes for 
the real part. 



1. Shape functions 



The shape functions are written 



^ K (7K(s),t) = 2 e 



7K 



Tif\/l+Oo|t| p -1k \A+aol*l 

- e^ K 



{2.1 



l + a \t\ ^A + a \t\ 

with the property ^ K (j K (s),t = 0) = 1. We have here introduced for each ampli- 
tude 4 energy dependent parameters, ax, (3r, Jk, ^k, whose roles are explained 
below. 7x is dimensionless, while ax, Ik and (5k are in GeV -2 . The fixed quantity 
ao = 1.39 GeV -2 is related to the square of the correlation length a of the gluon 
vacuum expectation value, with a = (0.2 ~ 0.3) fm [3]. 

These forms for the amplitudes originate in a formulation of pp scattering in impact 
parameter b space in studies based on the Stochastic Vacuum Model [H 0] . The 
analytic representation of our nuclear amplitude (except for the perturbative three- 
gluon exchange contribution) has a simple form in terms of impact parameter space 
through the Fourier transforms, 

T K (s,b) = ~ / d 2 qe-^ b T%(s,t = -g 2 ), 

Z7T J 

which are given in closed forms as 

T K (s,b) = -j-e ^ + \ K ij K (s,b) , (2.9) 
l p K 

where 

2 e~* K P ~\/ 7 ^ + ^o r _ Li , 62 

^ K ( Ss 6) = — , h2 e V 7 " +a ° (2.10) 

are the shape functions in 6 space. More details of the representation of the ampli- 
tudes in 6-space are given in Sec. |IV| We here choose to work with the amplitudes in 
t-space because they are more directly connected to the da/dt data that we intend 
to describe. 

These expressions are planned to represent the non - perturbative dynamics of scat- 
tering for all \t\ and the perturbative term R ggg is tuned to vanish for small \t\. 



The limits at \t\ = lead to the values for the total cross section a, the ratio p of 
the real to imaginary amplitudes, and the slopes Brj at t = through 



a(s) = 4^ (he) 2 + Aj(s)) , (2.11) 

Tg(s,t = Q) = a R (s) + X R (s) 
Tf(s,t = 0) aj(s) + Aj(s) 



P(-) = ^ "( = "T< T w\7 > (2-12) 



and 



dT#( 5 ,t) 



T#(s,t) dt t=o a K (s) + X K (s) 

'a K {s)(3 K {s) + ^A x (s)a (67x(s) + 7)] . (2.13) 

The determination of each amplitude starts with four energy dependent parameters 
that must be obtained from the data. However, with well established s dependences 
of the imaginary slope Bj (with typical Regge model behavior), connections between 
Xr+cxr and Xj+aj, plus constraints from dispersion relation for slopes [S], in practice 
we deal with fewer independent parameters to describe da/dt at each energy. 



As discussed before |3], the shape functions in Eq. (2.8) are very convenient choices 
for all values of \t\, determining consistently for all energies the zeros, the formation 
of dips and bumps, the signs and magnitudes of the two amplitudes and are able to 
reproduce with good accuracy all da/dt behavior. This description represents the 
non-perturbative QCD dynamics that is responsible for soft elastic hadronic scatter- 
ing. They account effectively for the terms of Regge and/or eikonal phenomenology 
that determine the process for \t\ ranges up to about \t\ ~ 2.5 GeV 2 . For higher \t\, 
perturbative contributions may become important. 



2. Universal behavior at large \t\ : Faissler measurements at 27.4 GeV 

For \t\ values beyond the dip and bump characteristic of the differential cross sections 
at ISR/CERN and Fermilab energies, namely for \t\ > 1.5 GeV 2 , the differential 
cross sections become increasingly independent of the energy. The measurements of 
pp scattering at y/s = 27 '.4 GeV [6] provides the only large \t\ data, covering the in- 
terval from 5.5 to 14.2 GeV 2 . It is known that the points are smoothly and naturally 
connected with the lower \t\ points at all energies [3]. As shown by Donnachie and 
Landshoff [7] this tail corresponds to a perturbative three-gluon exchange mecha- 
nism, with real amplitude contribution positive for pp and negative for pp scattering, 
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and da/dt falling as |t| -8 . The most remarkable example is given by the data of pp 
scattering at y/s = 52.8 GeV [5] , with measurements up to \t\ = 9.75 GeV 2 that 
superpose well with the 27.4 GeV tail. The observed behavior indicates that in this 
region the non-perturbative real part is indeed positive, and also that the magnitude 
of the imaginary part is small compared to the real part. At high energies in pp 
scattering the negative sign of the perturbative tail may lead to a marked dip in the 
transition region from 3 to 4 GeV 2 . 

Based on this expectation for the perturbative tail, we introduce a simple form of 
the universal three-gluon contribution amplitude R ggg (t) at large \t\, parameterizing 
the Landshoff term of the amplitude [7] as 

R ggg {t) = ±0.45 r 4 (l - e -° 005 l t l 4 )(l - e" ' 1 !*! 2 ) , (2.14) 

where the last two factors cut-off this term smoothly in the non-perturbative domain, 
and the signs ± refer to the pp and pp amplitudes respectively. 



Although the cut-off factors written in Eq. (2.14) have been adequate for all cases 
that were examined, their detailed forms in the transition range (2.5 < \t\ < 4) 
GeV 2 must be examined with data. 



B. Exponential forms of amplitudes for very small \t\ 

In the treatment of elastic pp and pp scattering in the forward direction, with 
amplitudes approximated by pure exponential forms, the differential cross section is 
written 

' h vr (he) 2 { e^m + p c {t) cos (a$) ] 2 



dt 



+ 



° e B lt /2 + pc^ sin ( a $)j 2 | ; (2.15) 



where t = —\t\ and we account for different values for the slopes Bj and Br of the 
imaginary and real amplitudes. The expression for the Coulomb phase that accounts 
for Br 7^ Bj derived in the Appendix. 

In elastic pp and pp scattering at all energies above \fs = 19 GeV, the real and 
imaginary amplitudes have zeros located in ranges \t\ ~ (0.1 ~ 0.3) GeV 2 and 
\t\ = (0.5 ~ 1.5) GeV 2 respectively, and the use of exponential forms beyond a 
limited forward range could lead to inaccurate determination of the characteristic 



FIG. 1: Analytic representation for the data on da/dt for elastic pp scattering at 7 TeV, 
together with the experimental data [2]; triangles are for subset A and open squares for 



subset B. The values of the parameters are given in Eqs. (3.1), (3.2). The inset at the 



up-right corner shows in closeup a gap between the datasets A and B, that seems to be 
at the limit of the error bars. 




forward scattering parameters a, p, Bj and Br. We then need to use more general 
expressions that are connected with the exponential behavior as limits, such as the 



forms with shape functions written in Eqs. (2.7) and (2.8). 



III. DATA ON pp SCATTERING AT 7 TEV 
A. Determination of Amplitude Parameters 

The TOTEM Collaboration has published data of differential elastic cross section 
[2] in two separate tables, in the ranges 0.00515 < \t\ < 0.371 (referred to as dataset 
A ) and 0.377 < \t\ < 2.443 (referred to as dataset B ) GeV 2 , with 87 and 78 
points respectively. Rather large systematic errors are informed in the dataset B. 



The representation of these data with the analytical forms of Eqs. (2.7) and (2.8) is 
shown in Fig. [I] The averaged square deviation for the dataset A + B (165 points) is 
(x 2 ) = 0.3105. The values of the quantities that also enter in the forward scattering 



amplitudes treated with exponential functions as in Eq. (2.15) are 



a = 98.65 ± 0.26 mb , p = 0.141 ± 0.001 , 
Bj = 19.77 ± 0.03 GeV" 2 , B R = 30.2 ± 0.7 GeV" 2 , (3.1) 
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FIG. 2: Data for elastic pp scattering at 7 TeV ([2], 165 points) and 52.8 GeV ( 0, 97 



points). The lines represent the parameterizations with Eqs. (2.7) and (2.8), with average 
squared deviations (x 2 ) = 0.3105 and 0.8328 respectively. Numerical characteristic values 
and description of features of the data and of the amplitudes are given in the text and in 
Tables U and U . 




TABLE I: Amplitude parameters for the 52.8 GeV and 7 TeV data. 





N 


a 


P 


Bi 


Br 


a i 


Pi 




Pr 


(x 2 ) 


GeV 


points 


mb 




GeV- 2 


GeV- 2 


GeV- 2 


GeV- 2 


GeV- 2 


GeV- 2 




52.8 


97 


42.49 


0.078 


13.04 


19.07 


5.9561 


2.3477 


1.1307 


1.1436 


0.8328 


7000 


165 


98.65 


0.141 


19.77 


30.20 


13.730 


4.0826 


4.7525 


1.4851 


0.3105 



and the other independent quantities that enter in the full-t forms are (all in GeV 

) 



cti = 13.730 ± 0.030 , /?/ = 4.0826 ± 0.0093 , 
\ R = 4.7525 ± 0.0155 , (3 R = 1.4851 ± 0.0318 . 



(3.2) 



Here we have 8 independent values, that are collected in Table [T] [23] Other quantities 



are related to these through Eqs. (2.11, 2.12, 2.13). Plots and features of the 



amplitudes are given in the next section. 
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Characteristic values of the differential and integrated cross sections for the 7 TeV 
data given by this representation are 

cr el . = 25.5418 mb ; a ine i. = 73.1082 mb ; o cX Ja = 0.2589 ; 

|t|di P = 0.4847 GeV 2 ; (da/dt) dip = 0.01532 mb GeV -2 ; 
Nbump = 0.6488 GeV 2 ; (da/dt) hmap = 0.02816 mb GeV" 2 ; 

ratio (bump/dip) = 1.8383. (3.3) 

In a closeup examination of the graph, shown in the inset in Fig. [TJ we observe a 
discontinuity in the magnitudes of da/dt in the junction of the two sets of data. We 
estimate that a renormalization factor 0.86 would adjust the higher \t\ to the more 
forward data in this region. 

It is very interesting to compare the 7 TeV data with the similar behavior of the 52.8 
GeV data [8], that are available in the ranges 0.107x 10" 2 < \t\ < 0.5546x 10~ 2 GeV 2 
with 34 points and 0.825 < \t\ < 9.75 GeV 2 with 63 points. This comparison is 
shown in Fig. [2] and in Tables [I] and [TT| and the behavior of the amplitudes is 
discussed in Sec. [TVl 



IV. AMPLITUDES 



The amplitudes obtained in the analysis of the data, based on Eqs. (2.7, 2.8), are 
shown in Fig. [3] and numerical information is given in Tables [I] and [TT]. Their general 
features are common to the lower energies from ISR and the Tevatron, with regular 
variation of the parameters. The results agree with requirements from dispersion 
relations for amplitudes and for slopes [5], and the Coulomb interference accounts 
for the necessary generalization of the Coulomb phase, presented in the Appendix. 
Near t ~ 0, the real part obeys the theorem by A. Martin [0] about its first zero, 
decreasing quickly and crossing zero at small |t| , before the imaginary part becomes 
small. The \t\ dependence of the amplitudes for all \t\ is shown in part (b) of Fig. 
[3j The imaginary part starts dominant over the real part, crosses zero at higher 
then remains negative and asymptotically tends to zero from the negative side, 
while the real part crosses zero again near \t\ = 0.8 GeV 2 , becoming positive. After 
the second real zero, namely for \t\ larger than about 1.0 GeV 2 , the real amplitude 
stays positive, without further oscillation, and for \t\ > 2 GeV 2 becomes increasingly 
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FIG. 3: (a) Forward scattering amplitudes Tr and Xj at y/s = 7 TeV in log scale, 
normalized to one at |t| = 0, showing their slopes, Bj = 19.77 and Br = 30.2 GeV~ 2 , 
and their curvatures, and indicating positions of the first zeros Zr(1) and Zi\ (b) Large-t 
dependence of the real and imaginary scattering amplitudes showing the complete set of 
zeros; (c) t dependence of the real and imaginary scattering amplitudes at ^fs = 52.8 GeV. 
Comparing with the figure for 7 TeV, we observe that all zeros move towards smaller \t\ 
values as the energy increases. 

dominant over the imaginary part. Important qualitative feature of our description 
of the data is that for large \t\ the magnitude of the imaginary part is smaller than 
the positive real part. 

This behavior is a regular continuation of the results obtained for the ISR energies, 
and as an example we show in part (c) of Fig. |3]the amplitudes for 52.8 GeV. At this 
and higher energies jl] the Zj zero occurs between Zr(1) and Zr(2). The magnitude 
of Ti(s,t) is in general dominant over the real part between Zr(1) and Zr(2), so 
that the dip is located near Zj. As the energy varies, while the first zero of Tr(s, t) 
stays almost constant in the \t\ range from 0.15 to 0.3 GeV 2 , the positions of both 
Zi and Zr(2) move to the left. The relative proximity of real and imaginary zeros 
influences the shape and depth of the dip. At 540 and 1800 GeV the imaginary zero 
is distant from both real zeros, so that no dip is formed, and only an inflection is 
observed in the da/dt curve. 

Table [IT] gives the positions of the zeros of the amplitudes and the characteristic 
observable quantities in pp elastic scattering at 7 TeV and 52.8 GeV. 

The details of the dip-bump structure carry information on the scattering ampli- 
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TABLE II: Positions of zeros of the amplitudes, locations of the predicted dip and bump, 
and ratio characterizing the shape of this structure. For comparison that shows the reg- 
ularity, we give the same information for ^fs = 52.8 GeV. We observe that as the energy 
increases all zeros move towards smaller \t\, with the imaginary zero moving faster, be- 
coming more distant from the the second real zero Zr(2) and closer to the first one Zr(1). 
*fs Zi Z R (1) Z R (2) |t| dip |t| bp (da/dt) dip (da/dt) hp ratio <r el cr ine i a c \/a 
GeV GeV 2 GeV 2 GeV 2 GeV 2 GeV 2 mb/GeV 2 mb/GeV 2 bp/dip mb mb 
52.8 1.3083 0.2710 1.6157 1.3560 1.7947 1.9 x 10" 5 6.3 x 10~ 5 3.2805 7.4308 35.0591 0.1749 
7000 0.4671 0.1641 0.8235 0.4847 0.6488 0.0153 0.0282 1.8383 25.54 73.11 0.26 

tudes, and are particularly sensitive to their relative behavior near their zeros [3]. 
Thus, a marked dip appears when one of the amplitudes (real or imaginary) crosses 
the zero in the interval where the other amplitude stays nearly constant with mag- 
nitude small compared to the variation of the former. In a domain where one of the 
amplitudes is dominant, the zero of the amplitude with smaller magnitude does not 
affect the observed differential cross section, as happens in the region of the first 
real zero (fa 0.15-0.3 GeV 2 ). 

The parts da R /dt and da 1 jdt of the differential cross section due to the real and 
imaginary amplitudes are shown in Fig. |4j 

Differently from lower energies, the dip at \t\ fa 0.5 GeV 2 is more influenced by 
the proximity of the imaginary zero (fa 0.5GeV 2 ) and the first real zero Zr(1) (at 
~ 0.2 GeV 2 ). Due to the increasing proximity of the imaginary and real zeros at 
higher energies, the dip/bump structure becomes more marked, and this is expected 
to happen at 14 TeV. 

At 7 TeV the bump that follows the dip is formed by the Tj and the Tr amplitudes 
that there have similar magnitudes, with Tj becoming more strongly negative while 
Tr quickly becomes zero. This explains the ratio nearly 2 (actually 1.838) between 
bump and dip heights. After the bump the fall is faster because it is only determined 
by the dominant magnitude of the imaginary amplitude, that decreases fast in \t\. Tr 
is small positive, but of long range, surviving until it meets the expected perturbative 
tail after \t\ = 2 GeV 2 . 
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A. Forward scattering 



The above description basically characterize our results of the analysis of all TOTEM 
data at 7 TeV, but the forward scattering part deserves a more detailed analysis, 
because the data in this domain have higher precision and are crucial for the deter- 
mination of the total cross section. 

To study separately the forward data investigating the use of amplitudes with the 



exponential forms of Eq. (2.15), we select a set of the first 40 points in the interval 
0.00515 < \t\ < 0.0907 GeV 2 . This forward set is used to compare the full-t 
analytic form described above with forward scattering analyses using exponential 



forms. Our analysis with Eqs. (2.7) and (2.8), with the parameters given in Eqs. 



(3.1) and (3.2), describes the forward set of 40 points with an average squared 



deviation ( X 2 ) = 0.1062. 



On the other hand, with Eq. (2.15) and the forward scattering parameters suggested 



FIG. 4: Partial cross sections da 1 /dt and da R /dt as functions of \t\ as calculated with 
the analytic forms of Eqs. (2.7) and (2.8). The dip in the sum da/dt (at 0.485 GeV 2 ) is 



close to the zero of the imaginary part (at 0.467 GeV 2 ). 
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by the TOTEM Collaboration [2] 

a = 98.58 ± 2.23 mb, p = 0.141 ± 0.007 , 

Bi = 19.89 ± 0.27 GeV~ 2 , B R = 5j , (4.1) 

we find (x 2 ) = 0.1129 for this set of 40 points. 

As a curiosity, we tested the pure exponential forms for the set of 40 points, using 



the values of a , p , Bj , Br given in Eq. (3.1), remarking that they are obtained in 
the all-t analysis. We then obtain an average (x 2 ) = 0.1034, which may even appear 
as a better description (actually, we can only say not worse) than obtained with the 



numbers of Eq. (4.1 ). 



Of course there is no statistical significance in these differences of (x 2 ) values, and 
visually these representations are almost degenerated as shown in part (a) of Fig. 
|5j The importance of the construction of the full-t analytical forms is that they are 
physically more realistic, including curvatures, magnitudes, signs, zeros of the two 
amplitudes, reflecting important information on the collision dynamics. In particu- 
lar, one should remind that the usual treatment of data with Br = Bj is essentially 
wrong, although it is used in practice because p is small. In fact, the slope B 
measured directly in da/dt is related to the slopes of the amplitudes through 

Bi + p 2 B R , . 

B = L + P p2 . (4.2) 

Therefore, the expected difference between values is approximately Bj « B — 
0.1 GeV~ 2 . However, such a simplification hides important properties of the dy- 
namics. Furthermore, since da/dt is not a pure exponential form A e~ B '*', determi- 
nations of a made with this approximation carry an intrinsic error. The scattering 
amplitudes must have zeros, consequently the slopes of the real and imaginary ampli- 
tudes are t dependent, deviating from pure exponential forms. So, the determination 
with high precision of total cross-section, p and slopes cannot be free from model- 
dependence, and differences are sensitive unless the measurements reach very small 
\t\ . Such a limitation of the approximation is shown in part(b) of Fig. |5j where 
we exhibit the large-t and the exponential solutions in the presence of the 87 points 
of the first part of TOTEM measurements, with \t\ < 0.371 GeV 2 . A pure expo- 
nential form (dashed line) appears as adequate (visually, and within errors) within 
the range \t\ < 0.25, but in reality the real part has a zero and becomes negative. 
We here stress that the physical quantities a , p , Bj and Br are obtained as lim- 
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FIG. 5: (a): Data and analytic representation for elastic pp scattering at 7 TeV in the 
forward range. Solid and dashed curves represent respectively the large-t solution and the 



exponential forms of Eq. (2.15) with parameters of Eq. (4.1) given by the experimentalists; 
(b) : The same functions are drawn together with the 87 points of first part of TOTEM 
data, showing the separation of exponential and full-t forms as \t\ becomes large. 





iting values as \t\ — > from a complex amplitude, thus requiring a basic analytical 
form. Therefore, they are naturally model-dependent, particularly if the data do 
no get close to the origin. It is thus very meaningful that the numerical values in 



Eqs. (3.1) and (4.1) are compatible. Note that our analytic representation for the 
amplitudes covers the whole t region with the same accuracy, and contains more 
physical information than a mere exponential fit. Naturally, the distinction between 
the two proposed descriptions for da/dt, and consequently for the values of a (98.65 
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and 98.58 mb) would require high precision in the measurements. For example, for 
\t\ = 5 x GeV 2 the predicted values for da/dt are respectively 447.2 and 444.8 
mb/GeV 2 . 

B. Universality of the perturbative three gluon amplitude at 7 TeV 

At energies above y/s ~ 20 GeV, a perturbative term due to three-gluon exchange, 
energy independent and of magnitude 10~ 7 to 10~ n mb/GeV 2 , produces a tail of 
form |t|~ 8 that dominates the scattering cross sections at large \t\. This term that 
was observed in Faissler experiment at 27.4 GeV [5J, is real positive for pp and 
negative for pp. In pp scattering at 52.8 GeV, where data exist up to \t\ ~ 10 GeV 2 , 
the matching of the data and the tail is perfect and shows that the real amplitude 
at intermediate \t\ values, namely after the dip in the cross section, should actually 
be positive, in agreement with our description of the amplitudes. 

In our approach, the non-perturbative contributions vanish fast for large \t\, and we 
expect that for \t\ ^> 3 GeV 2 , the cross sections at all energies behave like 

da/dt « (fic) 2 [Rgg 9 (t)} 2 « 0.08 r s (mb/GeV 2 ) . (4.3) 

Fig. [6] shows part of the 52.8 GeV and 7 TeV data, together with the points of 
the 27.4 GeV tail at large \t\. The superposition of points observed starting at 
\t\ = 2 GeV 2 leads to the feeling that with extended data the universal tail would 
be observed at 7 TeV, similarly to the behavior seen at lower energies. The second 
part of the figure shows the expected smooth connection of the TOTEM data at 7 
TeV to the 27.4 GeV tail data points, according to the construction of the ampli- 
tudes described above. According to this construction based on the assumption of 
universality, the \t\ distribution is expected to reach the three-gluon-exchange tail 
in the range from 5 to 15 GeV 2 . Unfortunately the experimental data end a little 
before, but observing the figure we are lead to believe that the present data are 
pointing towards the tail. Since T#(t) and the tail are positive, there will be no 
marked structure caused by cancellation, as expected to occur in pp. At large \t\ 
the cross section will be of form l/|t| 8 , according to the perturbative three-gluon 
exchange tail. The superposition between Tr an the perturbative tail may pro- 
vide precious information about the non-perturbative regime. The real tail, with 
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FIG. 6: (a) Prediction for the observation of the large \t\ tail (triangles) at 7 TeV, based 
on the comparison with the data at 52.8 GeV. The solid line is our solution for 7 TeV data 
with the tail added in the real part amplitude. In part (b) the dashed line is our solution 
without the perturbative tail. 




known sign, magnitude and shape, provides a reference basis, analogous to what 
the Coulomb interaction does at small \t\, important for the identification of the 
amplitudes. 



C. Impact Parameter Form for the Amplitudes 

As mentioned before, our amplitudes given as functions of the momentum transfer 



t\ by Eqs. (2.7, 2.8) are written in impact parameter space through the Fourier 



17 



transforms 



f K (s,b) = !-[ d 2 qe-^ b T»(s,t = -q 2 

where K = I and K = R in Tk(s, b) refer to the real and imaginary parts . They 
are given in analytical forms by 



* Pk 



where 



is the shape function in b space. Use has been made of the integration formula [4"l[T0] 

l w VT^ vdv =VWW (iA) 

Fig. [7] shows the parts T^ xp \s,b) and T^ hape \s,b) corresponding respectively to 
the terms in exponential and shape function of the amplitudes in \t\ , and Fig. [8] 
shows their sums. The shape functions are zero at b = and all parts are positive 
for all b. The amplitude in 6-space can be put in the form 

f R (s, b) + if ^8, 6) = i [1 - e 2i<5(s ' b) ] (4.5) 

with real S(s, b). 

The Fourier transform R ggg (s,b) of the perturbative tail is shown if Fig. [9J Its 
behavior is of a form approximated like R ggg (s,b) ~ 0.0022 J (c b)/(l + (d b) 2 ) , 
with c « 1.863 and d 4.0 in GeV. 



D. Comparison with the amplitudes in the BSW model 

Since the behavior of the amplitudes gives physical information on the collision 
dynamics and its determination is model dependent, it is important to compare 
our results with other models that describe elastic pp and pp scattering processes 
[TTHT5] . Results on amplitudes that can be directly compared with ours are given 
by the model proposed by Bourrely, Soffer and Wu (hereafter referred to as BSW 
model) [11]. We find important similarities and differences which we discuss below. 

The BSW model leads to values for the imaginary and real slopes with Br > Bj (as 
expected), to a first real zero at small \t\ (as expected) and to a first imaginary zero 
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near the cross section dip (as expected) and to a second real zero at intermediate 
\t\. These are the first basic, but crucial qualitative agreements with our approach. 
There also exist several quantitative and qualitative differences that may reflect 
important physical information behind the observables. In Fig. 10 we compare 
details of cross section (part (a)) and of amplitudes (part(b)). 

In the very forward region, the two descriptions cross each other, leading to different 
values at the optical point (with total cross section 92.33 mb in BSW and 98.65 mb 
for us) ; the slopes are qualitatively similar, with stronger Br leading to a first real 
zero at low \t\ in both expected by Martin's theorem. 

FIG. 7: (a) and (b) : respectively real and imaginary parts of the amplitude in impact 
parameter space , showing separately the contributions from the exponential ( Gaussian 
form, solid) and shape functions (dashed) of the original amplitudes, and their sum (dot- 
ted). 
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In the intermediate domain (0.3 < \t\ < 2) that includes the dip-bump region, the 
BSW cross section overshoots the experimental data and our results by a factor of 
nearly 2, caused by excessive magnitudes of both real and imaginary amplitudes 
in this region. The location of the dip-bump structure is also slightly shifted to 
a larger \t\ and the BSW dip is broader than ours. For \t\ ~ 2.5 GeV 2 , the two 
representations approach each other, but this is rather a coincidence as we can see 
from the behavior of the amplitudes. 

After the first zeros, the BSW amplitudes have larger magnitudes, and the second 
real zero is shifted to larger \t\ compared to ours, and this has influence in the 
dip-bump shape. Near the first zero of the imaginary amplitude, the BSW real 
amplitude is about 70% larger than ours in magnitude. 

Characteristic quantities of the BSW results are a = 92.33 mb, Bf = 19.3 GeV -2 , 
Br = 25.8 GeV -2 , p = 0.126, to be compared to the values given in Table | In the 
range of interest their amplitudes have two real zeros, at 0.205 and 1.475 (a third 
real zero is farther away), and three imaginary zeros at 0.525, 2.705 and 6.265 ; in 
our case there are two (and only two) real zeros, both more to the to the left, and 
only one imaginary zero, also at smaller \t\. 

FIG. 8: Real and imaginary parts of the amplitude in impact parameter coordinate. 
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Although the quantitative differences are important, it is impressive the similarity 
in the forms of the amplitudes in the two calculations in the data range, up to 2.5 
GeV 2 . This is the region where we believe that the non-perturbative dynamics is 
dominant and where eikonal-based models are expected to be valid. 

Important qualitative differences appear in the behavior of amplitudes for large \t\, 
namely \t\ ^ 2GeV 2 , which is illustrated in the form of partial contributions to the 



cross section of the BSW model in the up-right corner inset of Fig, 10 -(b). While 
in our model (solid for real and dash-dotted for imaginary), the non-perturbative 
amplitudes fall off exponentially for large \t\ and the the real part is made dominant 
with the action of the perturbative three-gluon exchange term, in BSW the ampli- 
tudes fall rather slowly, and exibit further zeros. Contrary to our approach, BSW 
cross section is dominated by the imaginary part for large \t\. Due to the presence of 
the second zero, there appears a slight inflection at \t\ ~ 3 GeV 2 . Interestingly such 
slight inflection also appears in our model coincidentally by a completely different 
reason, due to the perturbative real amplitude summing to the real amplitude. It 
seems that y/s = 52 GeV data supports strongly the contribution of this three-gluon 
exchange amplitude that we include in our description. 

FIG. 9: Representation of the amplitude for the perturbative tail in b space. 
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FIG. 10: (a) Comparison of da/dt at 7 TeV in our description (solid line) and in BSW 
model [11] (dotted line). In the large \t\ range our non-perturbative amplitudes fall ex- 
ponentially (dashed line) and the universal (for all energies) real three-gluon tail enters 
to dominate the scattering. In the BSW case a second imaginary zero at |t| ~ 2.7 GeV 2 
causes no mark, because the real part is not small, and a third zero in the imaginary 
amplitude predicts a marked dip at \t\ ~ 6.3 GeV 2 . (b) Comparison of amplitudes. Solid 
lines show the real and imaginary amplitudes of our representation, while dashed lines are 
for BSW. In the inset at up-right corner, we show the large \t\ behavior of the amplitudes 
of the BSW model in terms of partial differential cross sections; two imaginary zeros occur 
in this range. This can be compared to Fig. [4j 



The above described distinctions between BSW and our model lead to significant 
differences in the expectations for 14 TeV. In our case, we use the smoothness in 
the energy dependence of the parameters and of the characteristic features of the 
amplitudes to obtain a prediction for the higher energy, that shows a pronounced 
dip-bump structure, which appears at a bit smaller \t\ . This is due to the zero of the 
imaginary part, pinched by the two zeros of real part which are both approaching 



the imaginary zero for large -y/s. In Fig, 11, we compare the two predictions, which 
are quite different in the behavior of dip-bump structure, as well as in the large \t\ 
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FIG. 11: Comparison of predictions for 14 TeV , where the dip/bump structure becomes 
more pronounced. 

behavior. 

V. SUMMARY AND DISCUSSION 

In this paper we present an analysis of the recent TOTEM data for the elastic 
scattering cross section at y/s = 7 TeV, extending method employed at lower (ISR 
and Fermilab) energies [3]. The central point of our approach is that we describe 
the differential cross section directly in terms of scattering amplitudes, expressed 
analytically for all observed \t\ domain, together with the contribution from the 
perturbative three-gluon exchange process at large \t\. Differently from a mere 
parameterization fit of the differential cross section, our approach offers a qualified 
interpretation of data in terms of amplitudes, intending to serve as a bridge between 
data and theoretical descriptions by dynamical models. 

For example, it is fundamental to recognize the difference of the slope parameters of 
real and imaginary amplitudes to describe the all \t\ domain, although, due to the 
small value of the p parameter the effect is not clearly visible while looking only at 
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the forward scattering data. Furthermore, our analytic representation can also be 
written explicitly in the impact parameter space, that may permit us to discuss the 
physical meaning of our parameters in terms of geometrical models. 

Our approach works perfectly well for y/s = 7 TeV , giving an excellent rep- 
resentation for the data, including the forward region, and covering the dip at 
\t\ ~ 0.5 GeV 2 and the subsequent bump. The analysis leads to a description 
of the dip/bump structure of the differential cross section in terms of positions of 
zeros and magnitudes of the real and imaginary amplitudes. 

At 7 TeV the zero of Tj and the dip nearly coincide, at about 0.5 GeV 2 , so that the 
height of the dip informs the magnitude of Tr in this region. In the bump region 
after the dip, the magnitudes of the real and imaginary parts become similar, so that 
the form of the subsequent bump structure depends very sensitively on the interplay 
of the two amplitudes. The maximum of the bump occurs where the magnitudes of 
the real and imaginary amplitudes are nearly equal. With this behavior, the ratio 
bump/dip = 1.84 is nearly 2 . 

The gap in energy of the experiments at 7 TeV and the energies accessed in the long 
history of pp and pp scattering is enormous. Nevertheless, elastic (soft) processes are 
very conservative, with smooth energy dependences, so that the experience, methods 
and knowledge acquired before can be transferred to this higher energy domain. The 
zeros of the amplitudes determined from the present analysis (one imaginary zero 
and two real zeros) are consistent as natural extension of the previous work for lower 
energies. Since the behavior of the free parameters is smooth, we can extend the 
present analysis to the future \fs = 14 TeV experiment. We foresee that there will 
be a clear dip at t ~ 0.45 GeV 2 with similar bump structure with ratio ~ 3 and 
a total cross section a ~ 108 mb. This description concerns the non-perturbative 
regime characteristic of elastic and diffract ive scattering (including multi-pomeron 
and reggeon exchanges). 

It also becomes clear that the assumption, introduced and discussed by Donnachie 
and Landshoff |7j on the universality of the three gluon exchange amplitude, ob- 
served at 27.4 GeV [6] and higher ISR energies, is still consistent with the 7 TeV data, 
that likely are pointing to the distant tail, as seen in Fig. [6j If this assumption is 
true, some very interesting consequences can be drawn. First, for \t\ ^> 3 GeV 2 , the 
scattering amplitude is dominated by this process and consequently the imaginary 
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part should be small, while the magnitude and sign of the real part is determined 



by Eq. (2.14). Although the present analysis already indicates the consistency of 
this assumption, the precise behavior of the amplitudes in the transition region from 
non-perturbative to perturbative dynamics may be affected sensitively by the data 
in this region. Therefore, it is extremely important to confirm this behavior at 
y/s = 7 TeV by extending TOTEM analysis to higher \t\ values, with consequences 
for the identification of the amplitudes in the range of transition. If da/dt is found 
to follow the tail without marked structure, we can conclude that in this region 
the real amplitude has positive sign and that the imaginary part has much smaller 
magnitude. It would be nice to reach more extended \t\ ranges and high precision 
in this new era of studies of high energy collisions. 

It is interesting to remark that in the pp case, the negative perturbative term may 
produce a third zero in the real amplitude causing a dip in the \t\ distribution. This 
could have been seen in Fermilab data at 540 and 1800 GeV if the measurements 
had reached high enough \t\ values. 

In Subsec. IV D| we compare several aspects of our work with the results of a 
specific model [TT] . showing important similarities and differences, and stressing 
that comparison of models in the amplitude level is essential for the progress in the 
understanding of the scattering process. 

To conclude, working with a specific analytical form for the amplitudes, we have 
produced a detailed and precise description of the data of elastic pp scattering at 7 
TeV. The knowledge of the individual amplitudes carries more physical information 
on the dynamics of the scattering processes, and this work is part of an effort 
to find a consistent description of the amplitudes covering regularly the data for 
all energies in the whole t - range |4j. The present analysis of the TOTEM data 
reproduces accurately the behavior of the observed cross section in the whole t- 
range, with consistent proposal for the determination of the amplitudes, confirming 
the expectations of similarity with lower energies. It is hoped that this description 
can be successfully extended to the future measurements at higher energies. 
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VI. APPENDIX: THE COULOMB PHASE 



Here we present an expression for the Coulomb interference phase appropriate for 
forward scattering amplitudes with Br ^ Bj [16J. 

The starting point is the expression for the phase obtained by West and Yennie [17] 



= (-/+) 



In 



F N (s,t') 



F N (s,t) 



(6.1) 



dtf 

l-4p2 \f — t\ 

where the signs (— /+) are applied to the choices pp/pp respectively. The quantity 
p is the proton momentum in center of mass system, and at high energies 4p 2 « s . 

For small \t\, assuming that F N (s, t') keeps the same form for large \t'\ (this approx- 
imation should not have practical importance for the results), we have 

F N (s,t') 
F N (s,t) 

(6.2) 



F£ r (s,0)e s fl t '/2 +i F N ( s ,0)e s l t '/2 
F^(s,0)e B R t / 2 +i F^(sfi)e B l^ 2 



_c_ B R (t'-t)/2 , 
c+i c+i 



Bi{t'-t)/2 



where 



pe 



(B R -Bj)t/2 



(6.3) 



The integrals that appear in the evaluation of Eq. (6.1 ) are reduced to the form [18J 

rO A¥ r i 
1(B) = [ 1 - e B (*'-')/ 2 , (6.4) 



that is solved in terms of exponential integrals [19] as 

1(B) = #i[f(V + 

+ ln[- 



Bt] 
2 J 



f]+ln[f [Ap 2 + t 
f 2 7 . 



(6.5) 



The real and imaginary parts of the phase are then written 



*r(*> *) = (-/+) 





1 ■ 






+ c 2 + n 



c 2 I(B R ) + I(Bj) 



and 



*/(M) = (-/+)- 



7(5,) - I(B R ) 



(6.6) 



(6.7) 



With a in mb and t in GeV 2 , the practical expression for dcr/dt in terms of the 
parameters cr, p , Bj and Br is 



da If, \2 



+ 
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At high energies and small \t\ we simplify 4p 2 + t — > s and then the functional form 
of 1(B) is written 

+ 1b(-^+2 7 . (6.9) 

For large s , the term E 1 (Bs/2) can be neglected, and the phase becomes insensitive 
to s. 

The usual expression from West and Yennie 

Bt 



a^wY = (-/+)« 



7 + In 



2 



(6.10) 



can be obtained from Eq. (6.6) with Br = Bi = B and using the low t behavior 



rr Bt, , Bt 

H- -y] ~7 + log-— • 

We recall that the value of the slope B usually taken from the experimental dcr/dt 



data is the average given by Eq. (4.2). 



We obtain [16] that the values of a$>i(s, t) are very small so that the imaginary part 
of the phase can be safely put equal do zero. 

The construction of the Coulomb phase has been studied in the eikonal formalism 
appropriate for the interference of Coulomb and nuclear interactions [TS1 EDI EI] with 
special attention given to the influence of the proton electromagnetic form factor. 
These treatments keep the assumption that the real part of the nuclear amplitude 
has the same slope as the imaginary part, and the results are not very different from 
the West Yennie formula. Some of these results have been tested against the data 
E2l. 



Fig. 12 shows the values of the Coulomb interference phase in the kinematical 
conditions of the 7 TeV LHC experiment. We show together the phase obtained 
with the expression 



$cah„ = -[7 + H-Bt/2) + In (1 + 8/(5 A 2 )) (6.11) 
+ (-4t/A 2 )ln(-4t/A 2 ) - 2t/A 2 ] , 

where the proton form factor has been used with exponential form , A 2 = 0.71 GeV 2 , 
and B is taken as equal to Bj. 
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FIG. 12: Values of the Coulomb interference phase evaluated in the kinematical conditions 
of pp elastic scattering at 7 TeV , with account made for the difference in the slopes of 
the real and imaginary amplitudes. For comparison we show the phase obtained with the 
basic formula [20] of Eq.KU) 
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